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' We discuss the applicability, within the Random Matrix Theory, of perturbative 

treatment of symmetry breaking to the experimental data on the flip symmetry 
ly-^ . breaking in quartz crystal. We found that the values of the parameter that measures 

this breaking are different for the spacing distribution as compared to those for 
the spectral rigidity. We consider both twofold and threefold symmetries. The 
latter was found to account better for the spectral rigidity than the former. Both 
cases, however, underestimate the experimental spectral rigidity at large L. This 
discrepancy can be resolved if an appropriate number of eigenfrequecies is considered 
^ , to be missing in the sample. Our findings are relevant to isospin violation study in 

nuclei. 
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^ The study of wave chaos using acoustic resonators [H,^ supplies an invaluable addi- 

O . tional test of Random Matrix Theory (RMT) 0, 0|- In a 1996 paper EUegaard et al.^, 
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studied the gradual breaking of the presumed twofold flip symmetry of a quartz crystal 
by removing an octant of a sphere of an increasing radius at one of the corners and 
analysing the statistics of the resulting acoustic eigenfrequencies. They found a gradual 
. evolution of the spacing distribution from that of two uncoupled Gaussian Othogonal 
Ensembles (2G0E) when the crystal is an uncut perfect rectangle, into a single GOE, 
when a large chunk of the crystal is removed from one of the corners of the rectangle. 
^ ' This constituted a complete breaking of the symmetry present in the crystal in the un- 
. cut situation. The spectral rigidity, measured by Dyson's As{L) was also measured in 
this reference. The 2 uncoupled GOE's was found to underestime by a great amount 
the large-L data. This was attributed to pseudointegrable trajectories that do not suffer 
from the symmetry breaking. This point was further analysed by [lo|. Using techniques 



developed by Pandey [llj, Leitner [12| treated the symmetry breaking problem with 



RMT-perturbation. He addressed only the spacing distribution. This work was further 



extended to the spectral rigidity in [13|. In all of the above treatment of the data of 
jsl, the assumption was made that the uncut crystal has a twofold flip symmetry and 
thus is describable by two uncoupled GOE's. The treatment of Leitner [12j is found to 
describe fairly well the NNL distribution, but fails for the spectral rigidity, in contrast to 
the exact numerical simulation using the Deformed Gaussian Orthogonal Ensemble [5| 



recently performed in 1^. In this paper we further analyse the perturbative treatment 



of symmetry breaking within RMT. We find that the data of [8 



can be accounted for 



with SGOE's which are gradually mixed till a IGOE limit is attained. We further find 
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that that if some levels were missing in the sample of eigenfrequecies whose statistics is 
analysed, the A^lL) can be very well accounted for even at large L without the need for 
pseudointegrable trajectories, whose calculation is difficult. 



Using appropriate perturbative methods Leitner [12| was able to find a formula for 
the nearest neighbor distributio (NND) which contains the symmetry braking term. He 
started basically with the formula for the nearest neigbhour spacing distribution for the 
superposition of m GOE's block matrices 

where, for the case of all block marices having the same dimension one has 
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Emis)=iE,i-) , (2) 



E^{x)= / {1-F{t))dt, (3) 



X 



Fit) = / P,{z)dz. (4) 







In the above Pi{z) is the normalized nearest neighbour spacing distribution of one block 
matrix. It is easy to find for Pm{s), the following 

Pm{s) = - [{E,{s/m)r^' P,{s/m) + (m - l)(Ei(s/m))"^~2(l - F{s/m)f] (5) 

= p^J.\s) + Pl?:\s) (6) 

If all the block matrices belong to the GOE, then one can use the Wigner form for Pi{z) 

Pi{z) = -ze-^' ^-z, (7) 
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thus 



Er{z)=erfc\^z^ ^ 1 - z. (9) 

where the large-^ limits of Eqs. ©-([9]) are also indicated above. It is now clear that 
the above expression for Pm{s), ([5]) and ([6]), contains a term Pm\s) with level repulsion, 
indicating short-range correlation among levels pertaining to the same block matrix and 
a second term Pm\s) with no level repulsion, implying short-range correlation among 
NND levels pertaining to different blocks. Notice that for very small spacing, Pm{s) 
behaves as 

/ N TT m — 1 , , 

Pmis)-7^S + 10 

for m = 1, we get the usual -Pi(O) = 0, while for m > 1, we get Pm,(0) = (m — l)/m. 
To account for symmetry breaking, Leitner 13] considered the mixing between levels 



pertaining to nearest neigbhour block matrices and entails using the 2x2 P{s) distribution 
with full mixing. The DGOE result for the 2X2 matrix was derived in @ and the resulting 



P{s) is a product of a Poissonian term times a mixing term. Leitner's procedure 13 



amounts to multiply the factor Pm (s) of Eq. ([6]) by only the mixing term of the 2x2 



P(s) of [6|] with the mixing parameter A given by pjj], A = X^p^, with being the 
ratios of the variances of the matrix elements within a block matrix to that of matrix 
elements pertaining to neighbouring off diagonal block matrices., and p is the density of 
eigenfrequencies. Thus, he found, assuming that A << 1, 



PUs,A) = P^\s) + P2.2{s,A)P^\s) 



where P2x2(s, A) is given by [12| 

P2X2(S,A) = 
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where Jq is the modified Bessel function of order 0. Though Pm{s) is normalized, 
Pm{s,A) is not. Accordingly one supplies coefficients cat and cd, such that 



Pm{s, A, Cat, Cd) = CNPm{cDS, A) 



(13) 



is normalized to unity. Similarly, < s > should be unity too. Eq. (fTTj) can certainly 
be generalized to consider the effect of mixing of levels pertaining to next to nearest 
neighbour blocks, and accordingly, P3x3(s,A), given in Ref. would be used in Eq. 
( HT]) instead of P2x2('5,A). In the following, however, we use Eqs. (ITT]) . (IT^ as Leitner 
did pj- 



In Ref. [12|], Leitner also obtained approximate expression for the spectral rigidity 
A3(L) using results derived by French et al. Leitner's approximation to A3 is equal 
to the GOE spectral rigidity plus perturbative terms, that is 



Af)(L;A) ^ A3(L;oo) + 
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X ln(l 
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where 



TT 



2(r + 7r2A) 



(14) 



(15) 



For the cut off parameter we use the value [13[ r = CmG^' ^ , where c„ 
and 7 ~ 0.5772 is Euler's constant. This choice guarantees that when the symmetry is 
not broken, A = 0, A^"'^(L, 0) = mAsiL/m, 00). In Ref. [l^, Leitner fitted Eq. for 
m=2 to the NND from Ref. however, he did not fit the spectral rigidity. It is often 
the case that there are some missing levels in the statistical sample analysed. Such a 
situation was addressed recently by Bohigas and Pato [13] • These authors have started 
from the general expression of A3(L) derived by Dyson and Mehta namely. 
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where the two-point cluster function, Y2{xi,X2), which owing to translational invariance 
becomes a function of the difference a; =| — X2 |, is defined by the usual expression. 



where i?2 is the 2-point correlation function and R\ is the density of the spectrum. 

If a fraction, 1 — of the levels were actually analysed , the cluster function remains 
invariant, apart from a rescaling of the relevant variables, when the unfolded spactrum 
is employed, namely 

= 1 - n = ^2(^f , A\ (18) 

(1 -^)i?i(a;^)(l -^)i?i(x^) 
where the scaled variables x\ are just 



(1-9) 

Using the above equation for the cluster function in the general expression for A3(L), 



(19) 



we obtain the Missing-Level (ML) expression of [171 

A|(L)^4 + (l-9)'A,(^). 

In the application to our current problem of m-coupled GOE's, the above formula 
continue to be valid since the basic input into its derivation, namely the invariance of Y2, 
apart from the scaling of the argument x into x^, is quite general. Accordingly, we have 
the desired ML formula of A-^{L) for m-coupled GOE's, 



Ar^(L; A) = + (1 - g^Ar Aj . (20) 

The presence of the linear term, even if small, could explain the large L behavior of 
the measured A^lL). We call this effect the Missing Level (ML) effect. Another possible 
deviation of A3 from Eq. f[T^ could arise from the presence of pseudo-integrable effect 
(PI) 0, [l8|. This also modifies A3 by adding a Poisson term just like Eq. fl2U]) . 



The results of our analysis are shown in figures 1 and 2. In Fig. 1, the sequence of 
six measured NNDs were fitted for m = 2 and m = 3. It can be seen that the Leitner 
model with three coupled GOE's give a comparable and in some cases even better fit 
than the m = 2 one. Figure la in fact shows a rather sharp peak in our calculated P{s) 
for m = 3, P3(s, 0.0056). We consider this a failure of the Leitner formula (fT3!) for the 
uncut crystal. In fact, a more appropriate description of the uncut crystal is to take 
A = 0, namely a superposition of 3 uncoupled GOE's, which works almost as good as 
the 2 uncoupled GOE's description. The other parts of figure 1, (6) — (x) seem to show 
the same insensitivity of Pm{s,A) to m; the number of matrix blocks used in DGOE 
description. It is this insensitivity of the short-range nearest neighbour level correlation, 
measured by the spacing distribution, to the assumed symmetry inherent in the uncut 
crystal (and thus the number uncoupled GOE's employed to describe it) that forces us to 
examine the long-range level correlation, namely spectral rigidity, "measured" by Dyson's 
A3 statistics. 

In Fig. 2. the A-statistic was fitted with equation f|T4|) . It is clear from the figure that 
a good fit to the data of Ref. jl] is obtained with m = 3 for the values of A given in table 
1. This is to be contrasted with the case of m = 2 which, according to Eq. (HM results 
in A3^'*(L, A) that is always below the one with A = 0, A3^'*(L,0), which itself is always 
below the data points of Ref. [8|. For this reason, only the A3^^(L,0) is shown in the 



TABLE I: Values of A obtained by fitting Eqs. (1130 and (I14p respectively to the experimental 
NNDs and spectral rigidities from Ref. [g]. 
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A3(L) 


Data Set 


Ref. [15] 


Eq. (dlD m=2 


Eq. (dlD m=3 


Eq. m=3 


(a) 


0.0013 


0.0030 


0.0067 


0.0056 


(b) 


0.0054 


0.0063 


0.0098 


0.0016 


(c) 


0.0096 


0.010 


0.017 


0.0017 


(d) 


0.0313 


0.032 


0.064 


0.027 


(e) 


0.0720 


0.070 


0.13 


0.050 


(f) 


0.113 


0.12 


0.30 


0.16 


(x) 


0.138 


0.13 


0.34 


2.4 



figure. It should be noted that the A-statistics of the uncut crystal, Fig. 2a is very well 
described by that of 3 uncoupled GOE's, namely Af\L) = 3A3^''(L/3) which is always 
larger than the above mentioned A3^''(L) = 2A3^'*(L/2). The most conspicuous exception 
is Fig. [2]d which corresponds to r = 0.5 mm and where 1414 frequency eigenvalues were 
found. We consider this a potential ML case and take for A3, the expression given in 
Eq. (120|) and use it in Eq. (IT^ with g taken as a parameter. The results are shown in 
Fig. 3. We find perfect fit to the data, if g is taken to be 6%, namely only 94% of the 
eigenfrequencies were in fact taken into account in the statistical analysis. In contrast, if 
2G0E is used we still do not get very good agreement even if 18% of the levels are taken 
to be missing. There is, threfore, room to account much better for all cases (Fig. 2a, 
2c, . . . ) in the 3G0E description, by appropiately choosing the correponding value of g. 
This is not the case if a 2G0E description is employed. 

In conclusion, the perturbative tratment of symmetry breaking within RMT is assessed 
by using it to describe the coupling of m-fold symmetry The particular threefold case is 
used to analyse data on eigenfrequencies of elastomechanical vibration of a anisotropic 
quartz block. The treatment of Leitner [l2| is found to describe fairly well the NNL dis- 
tribution, but fails for the spectral rigidity, in contrast to the exact numerical simulation 
using the Deformed Gaussian Orthogonal Ensemble recently performed in [l^. However, 
by properly taking into account the ML effect we have shown that the A3(L), become 
closer to the data. 

We have also verified that if a 2G0E description is used, namely, m = 2 , then an 
account of the large-L behaviour of A3 can also be obtained if a much larger number 
of levels were missing in the sample. In our particular case of Fig. 2b, we obtained 
g = 0.18. This is 3 times larger than the ML needed in the 3G0E description. One 
major issue in the perturbative approach is the need to use different sets of values of the 
mixing parameter A for the NNL distribution, P{s), and for the A3(L). Further study of 
this approach is certainly required. Finally, we mention that the perturbation approach 
to symmetry violation study within RMT of the type discussed in this paper may be 



valuable to isospin breaking in nuclei. In Ref. [19l] the breaking of isospin was studied 
in the case of the nucleus ^^Al where both T = and 1 states are present in the low 
lying spectrum, which required the use of 2G0E description. Other, heavier, odd-odd 
nuclei may exhibit a spectrum where T = 2 states may also be present, requiring a 3G0E 
description of the isospin symmetry breaking. Work along these lines is in progress. 
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FIG. 1: Nearest Neighbour Distributions. Histograms show data (a)-(x) from Ref. j^. Thin 
and thick sohd hues show fits to the data carried out using Eq. (jl3p with m=2 and m=3 
respectively. In graph (x) the long-dashed line is the Poisson distribution, the dot-dashed line 
is the Wigner distribution and the dashed and dotted lines are the respective distributions for 
superpositions of 2 and 3 uncoupled GOEs. See Table I for the values of A obtained from the 
fits and the text for details. 




FIG. 2: Spectral Rigidities. Circles show data (a)-(x) from Ref. [8|. Thick solid lines show 
fits to the data carried out using Eq. (jl3p with m=3. The dot-dashed line is the GOE spectral 
rigidity and the dashed and dotted lines are the respective rigidities for superpositions of 2 and 
3 uncoupled GOEs. See Table I for the values of A obtained from the fits. 




FIG. 3: Missing Levels. Full line corresponds to A3(L) for the 3 GOE's of Fig 2b with g = 
0.06, while the dashed line corresponds to to the 2 GOE's with g = 0.18. The data points 
corresponding to the case with r = 0.5 mm, are from fig 2b of 0] 



